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Abstract. New spherically symmetric gravastar solutions, stable to radial 
perturbations, are found by utilising the construction of Visser and Wiltshire. The 
solutions possess an anti-de Sitter or de Sitter interior and a Schwarzschild-(anti)- 
de Sitter or Reissner-Nordstrom exterior. We find a wide range of parameters which 
allow stable gravastar solutions, and present the different qualitative behaviours of the 
equation of state for these parameters. 
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1. Introduction 

Black holes are widely accepted as physical objects due to their mathematical elegance 
and the strong astronomical evidence for their existence. The study of the properties 
of black hole horizons is of fundamental importance as horizons appear to provide a 
strong link between gravitation, thermodynamics and quantum theory. However, these 
horizons when described semi-classically give rise to a number of seemingly paradoxical 
theoretical problems which have yet to be satisfactorily resolved (see pQ for a review). 
For example, one such topic of current interest is whether a pure quantum state which 
passes over the event horizon of a black hole can evolve into a mixed state during black 
hole evaporation. This problem is known as the 'black hole information paradox'. It is 
generally accepted that the final resolution of this issue, and other difficult theoretical 
problems (such as the 'blue shift catastrophe' [2]) caused by the current description of 
black hole horizons, will be achieved using quantum gravity. However, we do not yet 
have a full theory of quantum gravity, and therefore, other possible solutions to remove 
the paradoxes generated by black hole horizons should be investigated. 

Given the above, it has been suggested that alternative endpoints of gravitational 
collapse of a massive star, which do not involve horizons, should be studied. The general 
idea is to prevent the possibility of a horizon forming, by stopping the collapse of matter 
at some radius greater than that of the horizon. This prevention of the formation of a 
horizon thereby precludes problems like the black hole information paradox. There 
are a variety of proposals including the Mazur-Mottola "gravastar" (gravitational 
vacuum star) |3j . Another scenario envisages the horizon as an emergent property of a 
Bose superfluid [4 . Mazur and Mottola's gravastar scenario is a solution of Einstein's 
equations which has a Schwarzschild exterior, a de Sitter interior corresponding to some 
sort of Bose condensate, and a rigid spherical shell of matter whose thickness is of the 
order of the Planck length, suspended approximately a couple of Planck lengths outside 
of the Schwarzschild radius. Due to their extreme compactness, it seems difficult to 
observationally distinguish gravastars from black holes. It has been argued that any 
star, with a surface, will emit much more radiation during accretion than black holes 
(which have no surface) (see [3] for a review), but it has also been shown that gravastars 
may be just as 'black' as black holes [§]. 

Although the gravastar model forwarded by Mazur and Mottola has not gained 
much attention by the majority of general relativists, it has led others to construct 
similar models. Bilic et al. consider a gravastar with a Born-Infeld-phantom interior 
geometry [2j , while Cattoen et al. generate a method for creating generalised gravastars 
with anisotropic continuous pressures jHj. 

Visser and Wiltshire [S| sought to determine whether the Mazur-Mottola gravastar 
is dynamically stable against radial perturbations. To do so, rather than considering 
the original Mazur-Mottola gravastar, they considered a simplified model with three 
layers: 



Stable gravastars with generalised exteriors 



3 



• An external Schwarzschild vacuum, p = = P 

• A single thin shell JU] 5 with surface density a and surface tension 9; with radius 
a > 2M. 

• A de Sitter (dS) interior, P = —p. 

Using this model they found a condition for stability for the thin shell against radial 
perturbations in terms of an 'effective energy equation' for a non-relativistic particle 
with 'energy' E in a 'potential' V(a), 

1 -d 2 + V(a) = E. (1) 

See section 4 from [5| for a precise definition of V(a) and the quantities therein. The 
thin shell will be stable against radial perturbations if and only if there exists some ao 
such that V(a) satisfies 

V(a ) = 0; V'(a ) = 0; V"{a ) > (2) 

Using Visser and Wiltshire's method one is able to prescribe the interior matter, m_(a), 
exterior matter, m + (a), and the potential V(a) of the gravastar to parametrically find 
the equation of state for the thin shell. Indeed, as we note below, the mass of the thin 
shell, m s (a) = 47r<j(a) a 2 , can be calculated as an explicit function of m + (a), m_(a) and 
V(a). 

Visser and Wiltshire demonstrated that there exist large classes of potentials, and 
consequently, equations of state, for which gravastars are stable against spherically 
symmetric gravitational perturbations as well as large classes of potentials which are 
unstable. Consequently, particular choices of potentials and mass functions need to be 
studied on a case by case basis, since more general criteria for stability are not presently 
known. 

In this paper we take the method as outlined by Visser and Wiltshire and analyse 
the stability of gravastars that have a Schwarzschild- (anti)-de Sitter or Reissner- 
Norstrom exterior. Given our current understanding, we pick these two types of exterior 
metrics as they are physically reasonable, spherically symmetric and static. 

2. Schwarzschild— (A) dS Gravastar 

2.1. Definitions 

Following directly on from the work of Visser and Wiltshire j^j, we start with the general 
equations, (47) and (49), from their paper. These relate the surface density, er(a), and 
surface tension, 9(a), of the thin shell which makes up the surface of the gravastar, to 
the potential and mass functions via 



(3) 



a(a) 



Ana 



1 - 2V(a) 



2m(o) 
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0(a) 



1 — 27(a) — m(a)/a — aV'(a) — m'(a) 
y/l-2V(a) -2m(a)/fl 



(4) 



where m_(a) is the interior mass profile, m + (a) is the exterior mass profile, V(a) is the 
effective energy potential, and 

[[X]] = X + - X. . (5) 

See ^T] for more details on this notation. See [5| for more discussion on the above 
expressions and definitions. We now specialise the exterior geometry to Schwarzschild- 
(A)dS, the interior to (A)dS, and the potential to zero. We therefore write 

m+(a) =M-Aa 3 /6, (6) 
m_(a) = ka 3 , (7) 
7(a) =0, (8) 

where A is the asymptotic constant spatial curvature of the exterior geometry and k is 
the curvature of the interior geometry due to a vacuum energy. We note that in this 
analysis, for E = in Eq.lfH). the choice of V(a) = provides a stable equilibrium 
for the gravastar thin shell, as a = 0, and hence the radius, a, cannot change to some 
other value. Simultaneously, V(a) = V'(a) = V" (a) = for all values of a. This 
specialization to V(a) = closely follows the calculation of Pj. We now convert all 
variables to dimensionless parameters, parameterised by M; 



A- — 

M 



L = AM 2 , K = kM 2 , p = aM, P = 6M . 



This leads to 
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(9) 



(10) 



Reality of the reparameterised surface density, p, and surface tension, P, requires 



2A 2 



6 3 



(12) 



2.2. Solutions 

We first present a detailed analysis of the dominant energy condition, as it provides 
the most restrictive conditions on the matter of the thin shell of the gravastar. Given 
that the matter is presumed to form through a quantum phase transition, we cannot 
specifically state the type of matter that forms the thin shell, but we shall assume the 
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dominant energy condition continues to hold. For brevity we define the region satisfying 
the dominant energy condition by Ai, and the boundary of Ai by dAi. The left and 
right sides of dAi are defined by the equations 

p + P = 0, p- P = 0, (13) 

respectively, satisfying the condition that p > 0. Solutions to the above equations give 
the points where the equation of state, p(P), enters or exits Ai. After some manipulation 
Eqs.(JT3J) become, 





2KA 2 = (3 - 8KA ) \ 1 - - + — - , (14) 



and 




(A - 3) VI - 2KA 2 = A\1-- + , (15) 



respectively. Squaring these two expressions can lead to an inclusion of a relative minus 
sign between the different sides of the equations, which we account for by placing extra 
restrictions on K and L. For K and L satisfying (fT2*|) we can see that the l.h.s. of (JUj) 
is positive. Therefore the r.h.s. is also positive, thereby providing the bound K < 
which applies to Eq.([16|) below. Similarly the r.h.s. of (fTB^) is positive, which implies the 
l.h.s. is positive and hence A > 3 for Eq. (jl7j) below. After some further manipulation 




Figure 1. The solution space for stable gravastars with (A)dS interior and 
Schwarzschild~(A)dS exterior satisfying p + P = in the left panel, or p — P = 
in the right panel. The number of solutions in each graph indicates how many times 
the equation of state intersects the part of dAi that the equation dehning the graph 
represents. The unshaded parts have no solutions, and correspondingly have no interval 
of p(P) which lies in Ai . 
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we find the multivariate polynomials, 

= (-32KL 2 - 192LK 2 ) A 8 + (lQL 2 - 576/T 2 ) A 6 
+ (240KL + 1152K 2 ) A 5 + (45L + 270K) A A 

+ (-324K - 120L) A 3 - 450KA 2 - 1084 + 225 , (16) 

and, 

= -124 + 27 - 6KA 4 + 3QKA 3 - 54/T4 2 - LA 4 . (17) 

The octic equation in A, (jlfijl . is not analytically solvable in general. We therefore solve 
this, and the quartic in A, (JT7J), numerically and create a contour plot representing the 
number of solutions for A in the phase space parameterised by K and L. Such solutions 
will have a number of restrictions placed on them. Specifically, they must satisfy the 
inequality (|T2"j) . p > 0, and any other conditions required for the consistency of (JTJj) 
and (JT3J). Figure H] displays the number of times the equation of state, p(-P), crosses the 
boundaries of M., p + P = and p — P = 0, separately. Figure |2 shows the combined 
results of the two panels in Fig. [T] Figure [T] is useful for determining the qualitative 
behavior of the equation of state, while Fig. |21 is useful for determining the existence of 
gravastar solutions which have a thin-shell satisfying the dominant energy condition. 




Figure 2. The solution space for stable Gravastars with (A)dS interior and (A)dS 
Schwarzschild exterior satisfying p + P = 0orp — P = 0. The unshaded region has 
no solutions where p{P) intersects dA4, and due to continuity of the the equation of 
state, no interval of p(P) lies in M. . The shaded regions with two solutions have one 
interval of p(P) which lies in Ai, while the region with four solutions has two intervals 
of p{P) satisfying the dominant energy condition. 

One can see from Fig. El that there are many different regions of K and L that 
have real and finite solutions to the equations that define dM.. This means that for 
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those values of K and L there is always an interval of the equation of state, p(P), which 
satisfies the dominant energy condition. This means there can always be gravastars 
with thin shells of matter which satisfy the dominant energy condition and are stable 
to radial perturbations. Visser and Wiltshire's gravastar with an exterior Schwarzschild 
solution is reproduced as a special case in this diagram when L = 0. The unshaded 
region has no solutions, implying that the equation of state does not enter or exit M.. 
One can check that for those parameter ranges p(P) does not lie in Ai. This means 
that for those values of the parameters K and L, it is not possible to form a gravastar 
which has a static thin-shell satisfying the dominant energy condition which is stable 
under radial perturbations. 
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Figure 3. The behaviour of p(P) for the thin shell of a gravastar with K = —0.01 
and a range of values of L. From left to right, L = {0.1, 0.06, 0, -0.0435, -0.7, -0.1}. 



One can find the five 'bounding curves', which are depicted by L(K) in Fig. |2j by 
studying the factorised discriminants of Eqs. ()lb| ) -([17| ) . where (|lb| ) -()17| ) are viewed as 
polynomials in A, and the restrictions on the parameters given by (|12j) that ensure p 
and P are real. The 'bounding curves' are the relevant parametric curves which bound 
the regions of the number of intersections of p(P) with dAi. They are given by, 

L= -1/9, 
L= — QK , 

L = 2K(V72K — 3) , (18) 

the ninth order polynomial in L, namely Lg(K), which is found by factorising the 
discriminant of (jlfij) . (The actual expression is too large to fit here.) 

There is a small region of the parameter space for K and L, with stable gravastar 
solutions when the interior is a de Sitter space, while there exists an infinitely large 
region of the parameter space for K and L when the interior is anti-de Sitter space. 



IK 
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The region with a Schwarzschild-de Sitter exterior is much larger than the region with 
a Schwarzschild-anti-de Sitter exterior. The most important bounds on K and L, 
governing the existence of gravastar solutions, are given by— 1/9 < L < L$(K). 

There are different qualitative behaviors that p(P) can take depending on the value 
of the parameters K and L, which we present in Figs. (0)-©. The various graphs show 
the quantitative behaviour of p(P) for specific values of K and L, but each graph is 
indicative of the qualitative behaviour of p{P) for a particular range of values of K 
and L. The equation of state exhibits a smooth transition from one type of qualitative 
behaviour to another as the values of K and L pass over the 'bounding curves'. 
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Figure 4. The behaviour of p(P) for the thin shell of a gravastar with K = and a 
range of values of L. From left to right, L = {0.07, 0, -0.0435, -0.07, -0.1}. 
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Figure 5. The behaviour of p(P) for the thin shell of a gravastar with K = 0.01 and 
a range of values of L. From left to right, L = {0.07, 0, -0.0435, -0.07, -0.1}. 
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Figure 6. The solution space for stable Gravastars with (A)dS interior and Reissner- 
Nordstrom exterior, satisfying p + P = in the left panel, or p — P = in the right 
panel. The number of solutions describes how many times p(P) passes through the 
left or right side of dAi in the left or right panel respectively. If there are no solutions, 
then p{P) does not cross dA4, and no interval of p(P) satisfies the dominant energy 
condition. 



3. Reissner— Nordstrom Gravastar 

We now consider a gravastar with charge q such that the radial electric field is given by 
E r = In this case, 

q 2 

m_ = kA 3 . (19) 

Following the same procedure as used for the Schwarzschild (A)dS case we find the 
dimensionless equations that define the endpoints of the interval of p{P) in A4, located 
on the left (p + P = 0) and right (p — P = 0) sides of dAi, given respectively by, 



3 " | + §) Vl^lKA 2 = (3 - SKA 2 ) + (20) 



and 



^i + ^vr^^-i + f, (21 ) 

where Q = > 0. The other dimensionless quantities are as shown in These can 
be manipulated to give 

= - QAK 2 A 8 + 128K 2 A 7 + (30K - QAQK 2 )A 6 - 36KA 5 - 20QA + AQ 2 
+ {2AKQ - 50K)A 4 + (40KQ - 12)A 3 + (3Q + 25 - 8KQ 2 )A 2 , (22) 
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and, 



- 2KA 6 + 12KA 5 + (-8KQ - 18K)A 4 + (2AKQ 
+ (3Q + 9 - 8KQ 2 )A 2 - Y1QA + 4Q 2 , 



A)A Z 



(23) 



respectively. Once again we demand that p and P be real which leads to 
1 



Q > 2A - A 2 . 



(24) 



As before, squaring (f2*Uj) and (j21j) can lead to an inclusion of a relative minus sign in 
(122)1 and ()2*3*)) respectively, which we account for by placing extra restrictions on K and 
Q so that that (|2*2j) and (J2SJ are consistent with (j20*)l and (J2TJ- We present our findings 
for the left and right sides of dAi, p + P = and p — P = 0, separately in the left and 
right panels of Fig. El while we combine the panels in Fig. Once again the parametric 



Two Solutions 




Figure 7. The solution space for stable Gravastars with (A )dS interior and Reissner- 
Nordstrom exterior, where p(P) for the thin shell satisfies either p+P = 0orp—P = 0. 
The regions with solutions have at least one interval of p(P) which lies in M. . Therefore, 
for these regions, there exist stable gravastar solutions. Conversely, the regions without 
solutions do not allow stable gravastar solutions where the thin shell satisfies the 
dominant energy condition. 

bounds are given by some of the factors of the discriminants of (122)1 and (|2H|) when they 
are viewed as polynomials in A, the condition that p > 0, and the inequalities in ([24)1 . 
The case studied by Visser and Wiltshire |JJ is recovered here with Q = 0. There are 
five bounding curves. Three of them are given by, 

Q = l, (25) 
3(2iT 2 ) 1 / 3 
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(26) 
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(27) 



The fourth one is given by polynomial which is twelfth order in Q, found by factorising 
the discriminant of ([20)1 . The fifth bound is given by the octic in Q found by factorising 
the discriminant of (|2ip. The qualitative behaviors of the equation of state depend on 
the choice of the parameters {K, Q}, as seen in Figs. IHlfTTI 
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Figure 8. The behaviour of p{P) for the thin shell of a gravastar with K = —0.5 and 
a range of values of Q. From left to right, Q — {0.9, 0.99, 1.01, 1.1}. The behaviour for 
Q < 1 is qualitatively different to that of Q > 1. 




Figure 9. The behaviour of p(P) for the thin shell of a gravastar with K = and 
a range of values ofQ. From left to right, Q = {0.99,0.999, 1.01, 1.1}. The behaviour 
for Q < 1 is qualitatively different to that of Q > 1. 



Extending this section to include a magnetic charge is straight forward, via the 
substitution q 2 — > q 2 + p 2 where p is the total magnetic charge of the system. It does 
not change the results in terms of Q, instead it amounts to a redefinition of what 
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Figure 10. The behaviour of p(P) for the thin shell of a gravastar with K = 0.1 and 
a range of values of Q. From left to right, Q = {0.9, 0.99, 1.01, 1.1}. The behaviour for 
Q < 1 is qualitatively different to that of Q > 1. 




Figure 11. This figure shows the behaviour of p(P) for the thin shell of a gravastar 
with K = 0.8 and a range of values of Q. From top to bottom, Q = {0.945,0.951, 1}. 
The qualitative behaviour of the equation of state is continuous over Q = 1. 



we attribute Q to. We note that including a non-zero charge to the exterior of the 
gravastar greatly increases the range of parameter K of the (A)dS interior for which 
solutions exist. Gravastars with a charge approximately equal to their mass are 'most 
favoured' for a de Sitter interior in the sense that for such values of Q the least amount of 
'fine-tuning' of K is required. However, when Q — 1, p(P) appears to be discontinuous 
for K < 0.423798525400 (13s.f.). 
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4. Discussion 

It has been demonstrated here that the method presented by Visser and Wiltshire 
can be used to generate stable gravastars with a non-trivial exterior involving either a 
vacuum energy or an electromagnetic charge. The Schwarzschild case they studied is 
retrieved here as a special case. The case with a Schwarzschild- (A) dS exterior, (see 
Fig. EJ), puts bounds on both the interior and exterior values of a vacuum energy 
allowable for stable gravastar solutions to exist. The most important bound in this 
case is that the exterior vacuum energy satisfies — \M~ 2 < A < L 9 (K)M~ 2 . 

The Reissner-Nordstrom solution (Fig. EJ) is very interesting as it is physically 
reasonable to allow massive stellar objects to have a (small) non-zero electric charge. 
We find that for a de Sitter interior, the range of the allowable vacuum energy smoothly 
increases as one adds charge, until the charge is close to the limit q 2 = M 2 . At this 
point the range of the parameter governing the interior vacuum energy becomes greatly 
enlarged by comparison with the q = (Schwarzschild) case. These results are important 
as they demonstrate that there exists a wide range of allowable gravastars with thin shells 
satisfying the dominant energy condition which are stable to radial perturbations. 

For future research it would be interesting to study the model in P with V(a) ^ or 
a ^ as these represent a wider class of configurations of potential gravastars. Another 
avenue for furthering this research would be to consider a gravastar with external Kerr- 
geometry, although the method provided by |Hj would have to be generalised first to 
include off-diagonal terms in the metric, and that is likely to be highly non-trivial. 

One source of a realistic gravastar interior might be a fundamental scalar field 
resulting from Kaluza-Klein compactification of extra dimensions ^2] (see or 
for a review) . Therefore it would be interesting to see what effect a scalar field has on the 
stability of the external geometry. Finally, investigating the implications of transitions 
in the equation of state from one type of qualitative behaviour to another is of interest 
(see Figs. I% HTU|) : we leave this work for a future paper. 
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